Chiral symmetry crossover with a linear confining 
and temperature dependent quark-antiquark potential. 
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Recently we developed a numerical technique to compute chiral symmetry breaking at T = 
with different current quark masses mo, including the current quark masses of the six standard 
flavours u, d, s, c, b, t. We also fitted from Lattice QCD data the quark-antiquark string tension 
a dependence on temperature T. We now utilize (j{T) to further upgrade the chiral invariant and 
confinement quark model to finite temperatures T 7^ 0. We study the quark mass at finite T and 
obtain the corresponding chiral crossover at T = Tc. The quark mass critical curve has a shape 
similar, but not identical, to the string tension critical curve. In the case of the lightest quarks, 
the quark mass and the chiral condensate essentially vanish at T = T^, except for the small explicit 
chiral symmetry breaking rriu and md- 



I. INTRODUCTION 



A. Motivation. 

We address the chiral crossover in the quark model 
perspective. The chiral crossover is a cornerstone to un- 
derstand the QCD phase diagram for finite T and 
fi. Notice that, after enormous numerical efforts, the 
analytic crossover nature of the finite-temperature QCD 
transition was finally determined by Y. Aoki et al. 0, 3 , 
utilizing Lattice QCD and physical quark, and reaching 
the continuum extrapolation with a finite volume analy- 
sis. The QCD phase diagram is extensively researched at 
the experimental collaborations LHC, RHIC and FAIR. 
This will help us to further understand how the universe 
evolves in the Big bang model of the cosmos. 

Here we utilize the linear confining potential for the 
quark-antiquark interaction, in the Coulomb Gauge Chi- 
ral Quark Model (CGvQM ), including both confinement 
and chiral symmetry [J]. While this model, inspired from 
the framework of Coulomb gauge Hamiltonian formal- 
ism is not yet full QCD, it is presently the only model 
able to include explicitly both the quark-antiquark con- 
fining potential and the quark-antiquark vacuum conden- 
sation. The CGxQM is for instance able to address ex- 
cited hadrons as in Fig. [T] and chiral symmetry at the 
same token, and it recently led us to suggest that the in- 
frared enhancement of the quark mass can be observed in 
the excited baryon spectrum at CBELSA and at JLAB 

Chiral symmetry breaking has been studied in detail 
with the CGxQM in the chiral limit and at vanishing 
temperature. It is quite well understood how, in the chi- 
ral limit of mo = 0, the quark develops a constituent 
running mass m{p) function of the momentum p. m(jj) 
is a solution of the mass gap equation (equivalent to 
the Schwinger- Dyson equation) for the quark, and this 
is important to understand the spontaneous breaking of 
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FIG. 1: First calculation of excited baryons with a chiral 
invariant quark model [^. 



chiral symmetry. Although the CGxQM is adequate to 
study the QCD phase diagram microscopically, the sci- 
entific community is only starting to explore |7H14|. the 
CGxQM with a finite temperature T ^ and with a 
finite current quark mass toq 0. Notice that a finite 
quark mass is not only crucial for the study of the hadron 
spectra, it is also relevant and for the study the QCD 
phase diagram. In the phase diagram, a finite current 
quark mass mo affects the position of the critical point 
between the crossover at low chemical potential fi and 
the phase transition at higher fi. The present work, not 
only addresses the QCD phase diagram, but it also con- 
stitutes the first step to allow us in the future to extend 
the computation of any hadron spectrum, say the Fig. [1] 
computed in reference Q, to finite T . 

We now review recent advances, leading to a temper- 
ature dependent string tension, and to a more efficient 
technique to solve the mass gap equation. These ad- 
vances are applied in Section II, where we derive the 
mass gap equation at finite Temperature. In Section III 
we solve the finite T mass gap equation, and compute 
the running mass m^ (p) to study the chiral crossover. In 
Section IV we discuss our results and conclude. 
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FIG. 2: We compare the T — static quark-antiqua rk p otential (solid line, left) to the T < Tc (left) and the T > Tc (centre) 
Lattice QCD data for the free energy _Fi, thanks to jl7H2ll |01af Kaczmarek et al. We also compare (right) the magnetization 
curve (solid line, right) with the SU(3) string tension critical curve (bullets) computed in Ref. [T^. We show (dashed line, 
right) the fit 1.21(1 - 0.990(r/rc)^)^''^ used in Ref. ^ to measure the finite string tension at Tc as an evidence for a first 
order phase transition. 



B. The quark-antiquark potential at finite 
temperature T 7^ 0. 

The most fundamental information for the quark- 
antiquark interaction in QCD comes from the Wilson 
loop and from the Polyakov loops in Lattice QCD, pro- 
viding the confining potential for a static quark-antiquark 
pair. This potential is consistent with the funnel po- 
tential, also utilized in the quark model to describe 
phenomenologically the quark-antiquark sector of meson 
spectrum, in particular to describe the linear behaviour 
of mesonic Regge trajectories. Notice that the short 
range Coulomb potential could also be included in the 
interaction, but here we ignore it since it only affects 
the quark mass through ultraviolet renormalization [l5| , 
which is assumed to be already included in the current 
quark mass. Here we specialize in computing different 
aspects of chiral symmetry breaking produced by linear 
confinement V ~ o'(T) r. 

The finite temperature string tension was computed 
for the first time by the Bielefeld Group for quenched 
SU(3) Lattice QCD [l^. The finite temperature static 
quark-antiquark free and internal energies at some finite 
temperatures T have also been computed in dynamical 
SU(3) Lattice QCD, by the Bielefeld for two light flavours 
Nf = 2 [IS-im, as depicted in Fig. H 

Recently we have empirically shown the string tension 
at finite T 7^ to be well fitted [l^, by a critical curve 
similar to the spontaneous magnetization of a ferromag- 
net [23} , i. e. solution of the algebraic equation. 
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dashed line at the right of Figure ^ , the difference be- 
tween these three scenarios is minute, and the validity 
of our fit is illustrated in Figure [21 Since the difference 
of these three scenarios has little effect in the numerics 
of the finite T mass gap equation, the empirical curve of 
Eq. (m provides us with the necessary function <t{T) to 
address confinement at finite temperature. 



C. The quark-antiquark potential, the mass gap 
equation and the Salpeter equation at T = 

To address the light quark sector it is not sufficient to 
know the static quark-antiquark potential, we also need 
to know what Dirac vertex to use in the quark-antiquark- 
interaction. This vertex is necessary to study not only 
the meson spectrum but also the dynamical spontaneous 
breaking of chiral symmetry. To determine what ver- 
tex to use, we review how the quark-antiquark potential 
can be approximately derived from QCD, in two different 
gauges. In Coulomb gauge 



V • A(x,i) = 



(2) 



T a{0) 



the interaction potential, as derived by Szczepaniak and 
Swanson [25l . [26| , is a density-density interaction, with 
Dirac structure 7*^(8)7". Another approximate path from 
QCD considers the modified coordinate gauge of Balit- 
sky and in the interaction potential for the quark sector, 
retains the first cumulant order, of two gluons (28l - [30j . 
This again results in a simple density-density effective 
7° (g) 7° confining interaction. Assuming such a Dirac 
^-^^ structure, our interaction potential for the quark sector 
is. 



Although the empirical curve of Eq. ([T]) corresponds to 
a second order transition (the solid line at the right of 
Figured]), and it is not a crossover as in full QCD, or a 
first order transition as in pure gauge SU(3) QCD (the 
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FIG. 3: (left) We plot the solution of the mass gap equation m(0) at T = 0, for different values of the current quark mass 
mo in the case where the string tension is cr = 1, and we also show, with a solid line, the fit with a two-parameter irrational 
function, (centre) We show a log log plot of minus the regularized quark condensate —{ijiil)) + (V''(/')o as a function of the current 
quark mass mo, with vertical dot-dashed lines representing the current masses of the quarks u, d, s, c, b, t. (right) We represent 
the regularized quark dispersion relation E{p) — ^ with increasing number of dashes per curve for five different current quark 

masses mo with respective values 10^^, 10~^, 10"^^'^, lO", 10^"'^. The quark condensate has an inflection point for finite quark 
masses close to the strange quark mass, and for large masses it rises linearly with mo- All results are in dimensionless units of 
a = 0.19 GeV^ = 1. 



where the Gell-Mann matrices are denoted A", and the 
density-density interaction includes just the linear confin- 
ing potential together with an infrared constant, which 
may be possibly divergent. While the CGxQM of Eq. 

is not exactly equivalent to QCD, we use it as our 
framework since it includes three crucial aspects of non- 
pertubative QCD, a chiral invariant quark-antiquark in- 
teraction, the cancellation of infrared divergences Isi]- 
Iscj , and a quark-antiquark linear potential p5l . [37l - l4l| . 
The mass gap equation and the energy of a quark are 
determined from the Schwinger-Dyson equation at one 
loop order using the hamiltonian of Eq. ([3]) . for a recent 
derivation with all details see 

In the limit of vanishing temperature T = 0, the inter- 
action in the four momentum of the potential and quark 
propagator term includes an integral in the energy 
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which factorizes trivially from the vector p momentum 
integral. Using spherical coordinates, the angular inte- 
grals can be performed analytically and finally only an 
integral in the modulus of the momentum remains to 
be computed numerically. We arrive at the mass gap 
equation in two equivalent forms, of a non-linear integral 
functional equation. 
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and of a minimum equation of the energy density £, 
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where the functions Ib and I a are angular integrals of 
the Fourier transform of the potential. In what concerns 
the one quark energy we get, 

(7 dk 

E{j>) = pC{p) + moS{p) + — / —lAik,p, m) X 

P Jo 

S{k)S{p) + lB{k,p,fi)C{p)Cik) . (7) 
We now extend this framework to finite T. 

II. IMPLEMENTING FINITE T 
A. The quark propagator at finite temperature 

r / 0. 

We now extend our equal-time density-density confine- 
ment to finite temperature. A framework for this exten- 
sion is the sum in Matsubara frequencies. With the equal 
time potential in our CGxQM , the integration in the 
three-momentum p and in the energy p*^ arc separable, 
and this is convenient to extend tour CGxQM with a 
Matsubara. 

We first study the wick rotation as a link between 
Eq. (g]) and the Matsubara sum. In the T = mass 
gap equation or Schwinger Dyson equation, we have the 
Minkowski space integral in p^ of the quark propagator 
pole of Eq. and this is equivalent to an integral in 
in Euclidian space after a Wick rotation in the Argand 
space. 
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This simple Wick rotation can be understood considering 
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FIG. 4: Artist's view of the Wick rotation. The Wick rotation 
from a real variable z = p° to a complex one z = ip^, can 
be performed without pole corrections, utilizing the depicted 
contour (solid in the figure), if there are only poles in the 
second and fourth Argand quadrants (dotted in the figure), 
and if the integrals in the circular paths of the first and third 
quadrants cancel in the limit of 1^1 — oo. 
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FIG. 5: The sum in Matsubara frequencies of the quark prop- 
agator as a function of -^j^ , including also energies below the 
fermi surface. The T^OoriJ — /.t— >oo limit of the sum is 
1/2, as in Eq. (g]) or in Eq. 



the closed contour of Fig. U 
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In the case there are no poles in the first or third quadrant 
and the circular paths cancel, one just has to replace p'^ 
by ip^ since in the real axis the path corresponds to 
z = p'^ and in the imaginary axis the path corresponds 
to z = ip**. 

Notice that the integral in Eq. is only identical to 
the one in Eq. ^ if the one quark energy E{p) > 0. If 
E{p) < the integral of Eq. ([5]) changes sign, and this is 
consistent with the pole moving from the fourth quadrant 
to the third quadrant of the Argand plane. While the 
integral in Eq. Q is insensitive to this translation of 
the pole, the contour of Fig. pleads to a pole correction. 
For simplicity, we choose to work with positive one quark 
energies only E{p) > 0. 

In finite temperature T and density p, the continuous 
euclidian space integration of eq. dS]) is extended to the 
sum in Matsubara frequencies, 
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i{2n+l)TTkT- [E{p) - fj] 



itanh^M^ 
2 2kT 



(10) 



where a chemical potential ^ may be included, as a finite 
density when a Fermi sphere of quarks is present. It 
is clear that in the vanishing temperature and density 
kT « E limit one gets back the initial Euclidean space 
integral of eq. ([5]) , when the Matsubara sum approaches 
the continuum integration with kT — >■ 41^. 

In what concerns the analytical calculations leading to 
the results in Eqs. (j4]) , ^ and (fTO)) . we used the even 
integral or summation to change the variable to a square, 
and we respectively used the residue theorem, real and 
indefinite integrals, and the analytical series summation. 



y 



■ tanh Try . 



(11) 



We verified all three analytical results with numerical 
integrations or sums. In Fig. ([5]) we plot the result of 
Eq. ([TU]) and it is clear that in the vanishing T limit all 
three analytical results coincide. 



B. Infrared regularization of the linear confining 
potential 



Notice that in the case of a linear potential, divergent 
in the infrared, the Fourier transform needs an infrared 
regulator A eventually vanishing. A possible regulariza- 
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FIG. 6: The linear potential must be regularized for the 
Fourier Transform, and we show our two different regular- 
izations of Eqs. (fT2|) and (fTl)) . both leading to V{r) — ^ ar 
plus a constant when r — >■ 0. Two successions of curves are 
plotted, to illustrate the limit of a vanishing infrared regulator 
A — >■ 0. The regularization of the negative curves maintains 
the potential monotonous and V{oo) = but adds an in- 
frared negative constant to the potential. The regularization 
of the positive curves maintains 1^(0) = but the potential 
decreases for large r. 



where the Fourier transform, 
V{k)=a -^'^ 



327rA2 
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(15) 



is such that the integrals in k no longer diverge. For 
instance. 



+00 



k^dkV{k) = 



(16) 



since this is proportional to V^(0) = 0. The new term 
in the potential ^^^^y^yi is equal to 5'^ {k) / \ in the 
limit A — > and thus the potential in Eq. ([T5|) is infrared 
finite. Both the potentials in Eqs. (|12p and arc il- 
lustrated in Fig. ini In the vanishing temperature limit 
T — Q the different regularizations lead to the same phys- 
ical results since any constant term in a density-density 
interaction has no effect in the quark running mass m{p) 
or in the hadron spectrum 

With the potential of Eq. (fT2|) the functions Ib and 
I A contributing to the mass gap equation and to the one 
quark energy are, 
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(17) 



tion of the linear potential is. 



C. A set of two non-linear equations. 
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corresponding to a model of confinement where the 
quark-antiquark system has an infinite binding energy 
J at the origin r = 0, is montonous and only vanishes at 
an arbitrarily large distance. This potential has a simple 
three-dimensional Fourier transform. 
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(13) 



and this is the most common form of the linear potential 
in momentum space utilized in the literature. Notice that 
this is infrared divergent due to the — 1/ A infinite bind- 
ing energy in the limit where the regulator A — >■ 0. If we 
want to avoid the infinite binding energy we should use a 
different regularization of the linear potential, also van- 
ishing when r — >■ 00 but not monotonous since ir grows 
linearly at the origin starting with 1^(0) — 0, 



V(r) 



ar e 



-Xr 



(14) 



At finite temperature T 7^ the mass gap equation for 
the quark running mass m(ji) couples to the one quark 
energy equation E{p) through the Matsubara sum, and 
thus we get a system of three non-linear coupled equa- 
tions. 
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lA{k,p, A) m^{k)p — lB{k,p, A) m^(p)k 
■\/A:2 + (fc)2 



p2 + mom^{p) a{T) 



^p2 ^ m^{pY 
lA{k,p, X)pk + lB{k,p, A) (p) m^(fc) 
^p2 jy^^^rpY^k"^ + m^{kY 

where the Matsubara sum normalized to 1 in the limit of 
small temperatures is, 



M^{p) = 2 ikT 

n—- 
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(19) 
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To solve the coupled systems of non-linear Equations 
(fT5)) . we apply the techniques detailed in Q, i. e. we 
utilize a Padc ansatz for the quark running mass (p) , 
we solve the mass gap variationally, we compute the one 
quark energy E^{p), we fit is with a Pade approximant, 
we feed it back into the mass gap equation, and then we 
repeat this cycle iteratively until the solution converges. 
The solution of the system of Eqs. ^TE\\ constitutes the 
main goal of this paper. 



III. THE T / SOLUTIONS OF THE MASS GAP 
EQUATION 

Finite temperature T changes the mass gap equation in 
two different ways, in the temperature dependence of the 
string tension cr(T) and in the Matsubara sum replacing 
the integral of the propagators. Since the finite tem- 
perature mass gap equation is numerically cumbersome 
to address, we first study separately the effect of each 
of these changes, before including both in the mass gap 
equation. 

Note that we work in units of ^/a ~ y/0.19 GeV, the 
string tension commonly used in the Charmonium spec- 
troscopy. For the transition temperature, we utilize the 
result computed in Lattice QCD di] of 0.176 (7) GeV. 
Thus we get Tc = OAO^. 
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FIG. 7: Here we consider only the effect of the temperature 
in the string tension (j{T), and plot the mass gap m(0) as a 
function of T for the light u (lightest) and d (slightly heavier) 
quarks. We have a crossover, close to a phase transition since 
the current masses of the light quarks are much smaller than 
the dynamically generated constituent quark mass. 
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A. The effect of the finite temperature string 
tension (t{T) 

Here we consider the T dependence of the string ten- 
sion only. Notice this approach is exact when our infrared 
regulator J7o — >■ oo as in the standard regularization of 
the linear potential of Eq. ([T^ . In that case E^^ (p) — > oo 
and the Matsubara sum A^^ = 1 is simply the same one 
occurring at T = . 

Then the finite T mass gap equation is similar to T = 
mass gap equation, just with a substitution of ctq by 
cr(T). The mass gap solution can be found with a simple 
rescaling of the string tension. 

In units of y/a = 1, we can fit the mass gap at the 
origin 77i(0) as a function of the running mass mp with 
a function interpolating from the chiral limit constant to 
the massive case linear function with a simple irrational 
ansatz. 



m(0) = - + ^'-o + Vi-^-mor_±^ (20) 

as depicted in Fig. ([3]- letf). Our fit produces a = 
-2.29223Vct, b = 0.656306v^. Then our solution is sim- 
ply found rescaling ^/a from 1 to ^/oiT) 

Our results, depicted in Figs. [7] and HI show a chi- 
ral crossover, since the mass gap never vanishes due to 
the finite current quark mass mg. We observe in our 
results that the heavier the quark, the weaker the chi- 
ral crossover gets, as show in Figs. [7] and [S] For the 
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FIG. 8: Log plot of the mass gap mq(0), considering T effects 
in (j{T) only, for the six different flavours of quarks as a func- 
tion of T. From bottom to top we show the quarks u, d, s, 
c, b and t. The heavier the quark, the weaker the crossover 
gets. 



heavier quarks the chiral crossover is very weak. On the 
other side, of light u and d quark masses, we still have 
a crossover since the mass gap starts finite at T = 
and continues to be finite beyond T^, although at large 
T >> Tc the mass gap is quite smaller than at vanishing 
T. 

Notice that the mass gap m(0) changes quite abruptly 
when T crosses Tc- The steepness of the light quark mass 
at T = Tc is due to our fit of (j{T) which is second-order 
like, as in Fig. [21 Possibly with a crossover-like <j{T) 
reflecting the dynamical fcrmion effects, the crossover in 
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FIG. 9: A possible scenario to compute the Matsubara sum 
consists in using the one quark dispersion relation subtracted 
by the enery at the origin E{p) — E{0) . This produces the 
smallest possible quark energy to be used in the Matsubara 
sum. The other opposite scenario consists in utilizing the 
energy E{p), wint an infrared divergence of In any case 
the one quark energy 
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the quark mass m(0) may get slightly smoother. However 
we do not know exactly what happens at the transition 
temperature for light quarks, since an opposite effect can 
make the string tension steeper again at Tc- Indeed for 
light quarks one may argue that the internal energy Ui 
(steeper than the free energy Fi) should be used. Thus 
in this paper wc cannot address in detail the exact tem- 
perature where T crosses but rather the remaining 
range of temperatures. 

B. The effect of the Matsubara sum effect in the 
propagator 

Here we are interested in finding the effect of the Mat- 
subara sum in the mass gap equation. To isolate this 
effect, we neglect any effect of the temperature T on the 
string tension, considering the string tension ctq sti all 
temperatures. 

We must then choose an infrared regulator, since 
adding a constant term —Uq to the potential, say — cr/A 
as in Eq. (fT2|) . affects the energy which is then shifted 
E^{p) E^{p) + Uo/2. Thus at finite T, the solution of 
the system of Eqs. (|18p depends on a constant shift of 
the potential, in contradistinction with vanishing temper- 
ature where M^{p) ~ 1 decouples the energy E^(p) from 
the mass [3| • At finite T the one quark energy (p) con- 



FIG. 10: Here we ignore the effect of the temperature in the 
string tension and consider only the effect of the Matsubara 
sum M^. We plot the running mass for the d quark , md{p), 
computed for the different temperatures present in Table [Jl 
ranging from T = to T = T^. Note that the temperature 
decreases the quark mass, but the effect is much weaker than 
the string tension effect in Fig. [T] 



tributes non- linearly to the Matsubara sum D^{p), and 
any change in i?^ (p) docs change the solution of the mass 
gap and energy coupled equations and thus the three Eqs. 
(dH]) are coupled. 

Here we choose a minimal Uq, closer to the regular- 
ization of the linear potential of Eq. ([T^ . just sufficient 
to cancel the one quark energy at vanishing momentum 
E^{Q) = 0. This case is also equivalent to the one when 
a chemical potential = E^ (0) is included, and thus 
E^{0) — fi = 0. We consider this case since it pro- 
duces the largest possible effect in the Matsubara sum 
Dt{p), interpolating between at vanishing momentum 
and 1 at large momentum. This does suppress the in- 
frared part of the integral if the mass gap equation up to 
E^{k)-E^{0) ~ kT. 

To exactly cancel the infrared divergences in the inte- 
grals, we utilize ansatze for the quark mass (p) and 
for the quark energy E^ (p) . Our numerical technique 
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FIG. 11: Full computation of the mass gap for the d quark 
, Tnd(O), in the infrared regulator case where the Matsubara 
sum effects are maximal. Although not as steep close to T = 
Tc, the curve is quite similar to the one where we ignore the 
Matsubara sum effects. Again we have a crossover at Tc. 



was tested in great detail for T = an here we utilize the 
same ansatze and numerical techniques Our ansatze 
for (p) and (p) are the Pade approximants, 
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■K p^ — ■m^(p)'^ 

+pC{p)+moS{p) + E^{p) , 
1 



Co + eip + e^p-^ + e^p^ 



After each iteration the resulting [p) and Er^ are fitted 
again with our Pade approximants, and a new iteration 
is started. But for tch solution for each temperature and 
mass is completely determined by the three parameters 
for the fits of m^{p), since E^{p) is a function of the 
quark mass. 

To ensure convergence, we utilize the technique Q of 
translating the mass gap equation into the variational 
equation for the vacuum energy density. 



£ = 



dp 
2^' 



2tt 



D^{p)\2p^C{p)+2p^moS{p) 



lA{k,p,fi) S{k)S{p) 



+lB{k,p,pi)C{p)C{k) 



(21) 



now with a new dependence on temperature in (mo- 
mentarily we consider ut = ctq) and in the Matsubara 
sum AA^ . When the current mass is small compared 
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FIG. 12; Full computation of the running mass for the d quark 
, md{p), computed for the different temperatures present in 
Table III ranging from T = to T = Tc. Unlike when the 
Matsubara sum only is considered, here the quark mass indeed 
equals the current quark mass Md, almost vanishing, for T = 



with the typical scale ^/a(T) of our problem, we uti- 
lize directly this equation. In the opposite case were the 
string tension is much smaller than the current mass, 
close to Tc, then the solution m{p) is close to the current 
mass toqi 



m{p) ~ 7710 -I ) , 

7770 mo 



(22) 



and then we only need to compute the dimcnsionless 
function J-{p), produced pcrtubatively Q after few it- 
erations of the fixed point equation. 

We consider the case of 777o = 777^ corresponding to a d 
quark. Our results are illustrated in Fig. (TUl Note that 
the effect of the Matsubara sum TW^ is to decrease the 
quark mass with increasing temperature, but neverthe- 
less the Matsubara sum effect in Fig. [10] is much weaker 
than the string tension effect shown in Fig. [T] 



9 



mo 



T 



0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 
0.00527656 



0.000000 
0.0313837 
0.0625738 
0.0933782 
0.123607 
0.153073 
0.181596 
0.208999 
0.235114 
0.259779 
0.282843 
0.304162 
0.323607 
0.341056 
0.356403 
0.369552 
0.380423 
0.388948 
0.395075 
0.398767 
0.400000 



2 



Co 



C2 



C4 



-0.558232 
-0.514877 
-0.482888 
-0.456189 
-0.425455 
-0.390516 
-0.35179 
-0.309921 
-0.265958 
-0.221304 
-0.177546 
-0.136292 
-0.0990367 
-0.0670518 
-0.0412995 
-0.0223216 
-0.0100832 
-0.00365887 
-0.00105524 
-0.000187953 
-0.0000000 



6.20966 
6.24603 
6.29268 
6.35722 
6.45135 
6.57955 
6.74966 
6.97245 
7.26188 
7.63714 
8.12658 
8.77490 
9.65716 
10.9094 
12.8039 
15.9619 
22.0570 
36.2780 
73.6460 
190.883 
oo 



25.2077 
25.1480 
25.0808 
25.0130 
25.0128 
25.2058 
25.6934 
26.5531 
27.8649 
29.7358 
32.3267 
35.8925 
40.8552 
47.9602 
58.6492 
76.1081 
108.769 
183.098 
382.799 
1073.65 
oo 



14.5339 
14.6870 
14.8787 
15.1468 
15.5903 
16.3132 
17.4403 
19.1262 
21.5887 
25.1689 
30.4365 
38.3987 
50.9525 
71.9645 
110.150 
187.827 
371.573 
886.265 
2491.37 
9021.24 
oo 



TABLE I: The parameters for the quark functions m^, (p) 
defined in Eq. (pTj) . 



C. Full solution of the mass gap equation at finite 
temperature T 

Here we determine the full solution of the mass gap 
equation at finite temperature T, utilizing Eq. (|2ip . but 
combining the temperature effects in the string tension 
<t(T) with the temperature effects in the Matsubara sum 
Mj.{p). In the present case we continue to consider that 
fi ~ E{0), the case maximizing the finite temperature 
effects of the Matsubara sum. 

We now fix the current quark mass to mg ~ md, quark 
d current mass and run the temperature T so that we get 
a curve. We numerically consider a set of temperatures 
denser at Tc, with the parametrization of T = Tc cos 
where the 6 considered are equally spaced in the interval 
ranging between tt/2 and 0. We show the best fitting 
parameters cq , C2 , C4 , in Table HI we depict the mass gap 



critical curve in Fig. [TT] and we show the different run- 
ning quark masses in Fig. [T^] . It occurs that indeed 
a crossover is found in the mass gap critical curve at 
T = Tc, and that the critical curve is similar (but not 
identical) to the critical curve computed with the string 
tension temperature effects only. 

IV. CONCLUSION 

We apply a new variational technique, in the frame- 
work of a CGxQM , to solve the mass gap equation at 
finite temperature T. The quark mass m^{p) and the 
quark energy E^^ (p) are fitted with Pad approximants, 
the quark mas parameters are displayed in Table H] 

Wc compare the two different temperature contribu- 
tions, in the string tension cr(T) and in the Matsubara 
sum M^{p). It occurs that the dominant contribution, 
and also the simplest one to apply in the quark model, 
is the one of the string tension cr{T). This happens 
not only because the dcconfinemcnt critical temperature 
Tc — QAQ^/a is relatively small when compared to the 
string tension at vanishing temperature (Tq, but also be- 
cause the Matsubara sum never vanishes while the string 
tension really vanishes when T — >■ T^. Thus the quark 
mass critical curve has a shape similar to the string ten- 
sion critical curve, but the curves are not exactly identical 
since the quark critical curve is less steep at T ~ T^. 

Moreover, since the light current quark masses rriu and 
rud are small compared with the string tension at vanish- 
ing temperature ctoj the quark mass gap m(0) essentially 
follows the string tension curve. This is a quite simple 
result, relevant for further studies at finite temperature. 

With our excellent fits of the dynamical quark mass 
rriy {p) and of the quark dispersion relation E^^ (p) we are 
now well equipped to address further studies, such as the 
hadronic excited spectra at finite temperature T. 
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